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0^ ' We consider a cosmology with a noncompact nonlinear sigma model. The target space is of de 

Sitter type and four scalar fields are introduced. The potential is absent but cosmological constant 
term A is added. One of the scalar fields is time dependent and the remaining three fields have no 
time dependence but only spatial dependence. We show that a very simple ansatz for the scalar 
O I fields results in the accelerating universe with an exponential expansion at late times. It is pointed 

D , out that the presence of the energy density and pressure coming from the spatial variation of the 

• three scalar fields plays an essential role in our analysis which includes A = as a special case 
and it discriminate from the standard A-dominated acceleration. We perform a stability analysis 
of the solutions and find that some solutions are classically stable and attractor. We also present a 
nonperturbative solution which asymptotically approaches an exponential acceleration and discuss 
possible cosmological implications in relation to dark energy. It turns out that the equation of state 
approaches asymptotically uj — — 1 both from above and below, but the crossing does not occur. It 
predicts present value of a; ~ — 1 T 0.07, which is within the region allowed by the observational 

»' | ' data. This solution also exhibits a power law expansion at early times, and the energy density of 

the scalar fields mimics that of the stiff matter. 

PACS numbers: ll.10.Lm, 95.36,+x, 98.80.-k 
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CN ' I. INTRODUCTION 

The recent cosmological observations [l| provide many precise data and arouse an explosion of recent interests in 
the cosmology. The most recent data and its cosmological interpretation [2j indicate that about 73% of our Universe 
is made of dark energy, the origin of which is one of the greatest puzzles in the modern cosmology [3j . 

It is highly conceivable that the dark energy is responsible for the late acceleration of the Universe [i[ and many 
candidates have been proposed. The simplest approach for the accelerating universe is to introduce the cosmological 
constant 1 1, 1 ,1 1 (1 for the dark energy. Other approaches [4[ include dynamical models of the cosmological constant 
[ToL ITU [la liallll. [Ta . El El EE El, ill EL [H, lM HI • Among them, most commonly proposed candidates are the 
quintessence, which is described by a scalar field minimally coupled to Einstein gravity with a potential [l2[ EB|- It 
is shown that the scalar energy density is subdominant in the matter dominated, and then, acceleration takes over 
at later stage of the cosmological evolution. Later the phantom model with a negative kinetic energy scalar field was 
proposed [2 a] to account for the region where the equation of state is less than u = — 1, and quintom model where the 
ordinary scalar and the phantom are both introduced 1261 to explain the crossing of the uj = — 1 line. In these models 
and subsequent works, various forms [H El El El El 13 of the potential are introduced to account for the dark 
energy and the asymptotic acceleration is achieved through the quintessence or phantom dominance at late stage. 

In this paper, we consider a cosmological model with a nonlinear sigma model [28| with a cosmological constant term. 
The target space is noncompact and is of four-dimensional de Sitter manifold and four scalar fields are introduced. 
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One of the motivations is that the scalar fields have a geometric origin and the potential term is not necessary. Another 
motivation is to consider the spatial dependence of the scalar fields and examine its consequences. To solve the Einstein 
equation, we assume that only one of the scalar fields is time dependent and the remaining three fields have no time 
dependence but only spatial dependence. We first show that a very simple ansatz can solve the Einstein equation 
which describes the acceleration of the universe with an exponential expansion at late stage. The spatial contribution 
in combination with the cosmological constant forms an effective cosmological constant and plays important roles in 
order to provide the necessary energy density and pressure. We will also show that the acceleration is possible even 
without the cosmological constant. It seems that this feature of contributing the energy density and pressure coming 
from the spatial variation of the scalar fields was not considered before in relation with dark energy. 

Then, we perform a stability analysis of the solutions. We find that some of the solutions, depending on the values 
of the parameter given, are classically stable and attractor solutions. We will consider two cases where the target 
space has signature (+, — , — , — ) or (— , +, +, +). Especially, in the (+, — , — , — ) case, the linear stability analysis fails, 
but we are able to find out that there exists a nonperturbative solution which asymptotically approaches the dc Sitter 
acceleration, but at early times it is a power law expansion. In this case, the cosmological constant term is uniquely 
fixed in terms of the other parameter. These features differentiate the present analysis from the standard A-dominatcd 
late-time exponential acceleration. 

One might think that adding a cosmological constant term with the scalar fields could be ad hoc, but it seems 
that at present, the dynamical models of the dark energy is not completely successful in solving the cosmological 
constant problem and many of them require some kind of fine tuning anyhow. Nevertheless, the aim of this paper is 
not to explain the smallness of the cosmological constant (the fine tuning problem in our approach is mentioned in the 
Conclusion and Discussion section) , but to focus on the late-time exponential acceleration of the universe and stability 
of its behavior. It turns out that the accelerating universe requires some bound on the original cosmological constant 
term. In (— , +, +, +) case, it require that the original cosmological constant term must be negative for stability, still 
the acceleration is possible and can be led by the scalar fields. 

The paper is organized as follows. In Sec. II, we present noncompact nonlinear sigma model coupled with Einstein 
gravity with a cosmological constant term and discuss the ansatz which solves the equations in some generality. In 
Sec. Ill, we describe the exponential accelerating solution with de Sitter target space. In Sec. IV, the stability 
analysis is performed and allowed range of the cosmological constant is classified. In Sec. V, a nonperturbative 
solution is obtained and possible cosmological implications in relation to dark energy is given. Section VI includes 
the conclusion and discussion. 



II. THE ACTION AND COSMOLOGICAL CONSTANT 



We consider an action in which the Einstein gravity is coupled to a nonlinear sigma model with a cosmological 
constant term (in units of M p = 1): 

S = J d 4 xV^9 [ \R - 9 -^G af3 {^)d^ a d v ^ - A + £ matter } (I) 

where $ Q = (0, a 1 ) (i = 1, 2, 3), G a p is the metric of the noncompact target space, A 2 is the self-coupling constant of 
the nonlinear sigma model and it is assumed to be positive. A is the cosmological constant. The equations of motion 
are given by 

Rnv = ^G a pd^ a d^ + Ag^ + % v (2) 
1 --d,[^—gg^G aP d^] = l^lg^d^d^, (3) 



—-ml* ^ j 2 9$Q 

where T^ v (= — g^T/2) is assumed to take the perfect fluid form; 

= (-pm, p m , p m , Pm)- (4) 

The matter sector satisfies the continuity equation; V ^T^ v = 0. 

If we ignore the matter part, we can solve the Eqs. @ and ([3]) with the following ansatz 

0= t, a 1 = x\ (5) 
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To check whether Eq. ([2]) can be solved (without T^) with this ansatz, first note that if A = 0, g^ v {t, x l ) = +G llu ((j), a 1 ) 
satisfies the equation as long as the scalar curvature of the space-time metric g^ v and that of the target space G a p 
are constants. Then, we can add an cosmological constant A which has the same sign with these curvatures and the 
equation can still be satisfied [2!|[3(|. This has the effect of scaling the space-time metric via g^ — ► (1 + 4A/i?)g /i „, 
where R is the scalar curvature constant. The cosmological constant could even have some value of the opposite sign 
with these curvatures as long as the absolute value of the cosmological constant is smaller than that of the target 
space, i.e., |A| < |i?|/4. This point can be extended further. Suppose g^(t,x l ) = — G l _ LV ((f) 1 a 1 ). Then, the scalar 
curvatures of the space-time and target space have opposite signature and the equation cannot be satisfied with A = 0. 
But if we add A such that the sign is the same as the space-time scalar curvature constant and the absolute value 
is greater than the scalar curvature of the target space, the equation can be satisfied. Also, one can check that the 
metric ansatz g^ v {x) — ± G^{4>,a l ) satisfies the Eq. ©. In summary, we find that the metric ansatz and © solve 
© and and we have 

R llv = (±-^ + k)g flv , (6) 

as long as the constant curvature condition is satisfied and without T^ v . 

The a 1 = x 1 ansatz first appeared in higher dimensional gravity theory in association with spontaneous compactifi- 
cation of the extra dimensions [29|, [3(| ■ It does not break the isotropy and homogeneity of the universe as long as we 
do not introduce the potential for the a fields. Also the (f> = t has been exploited to unify early-time and late-time 
universe based on phantom cosmology [3l|. Note that the quantity A e // = ±2/ A 2 + A plays the role of the effective 
cosmological constant and there is curvature constant restriction on the value of A; ±2/A 2 -I- A must have the same 
signature as that of the space-time scalar curvature constant. The above aspect of the ansatz (O, ([6]) is quite a general 
feature of the nonlinear sigma model coupled to gravity. In this paper, we will consider the de sitter target space with 

GS = e(l, -e 2 ^, -e 2 ^, - e 2 ^ ) (7) 

with e = =Fl and £ being an arbitrary positive constant. It turns out that the spatial ansatz ([5]) provides contribution 
of the energy density and pressure such as to reveal diverse aspects of the late time exponential acceleration, not 
present in the standard cosmological constant dominated acceleration. We will also find that the allowed value of the 
cosmological constant divides further if required the stability. 



III. DE SITTER SOLUTION 



To discuss the cosmological implication of the solution ([5]), ([6]) with the de Sitter target space metric ([7]), we 
introduce the standard space-time metric via 

ds 2 = -dt 2 + a 2 (t)dx l dx l . (8) 

With H — d/a, the equation of motion © becomes 

B = -^4 4- 2 -^*^-^{Di4') 2 -^e 2 f*(d i <TJf]-l(l+a, m )p m , (10) 

where u) m — p m /Pm- The continuity equation implies p m cx a - 3 ( 1 +"™). Plugging the ansatz a 1 = x 1 and <j) = <^(£) 
the above equations, Eqs. ©, © and pH|) become 

= + 3H0-3^— , (11) 
a z 

H = — h «[ - -W 2W 1 - + (13) 

The second terms in both (|12[) and (|13p are the contributions coming from the spatial variations of <ii which is essential 
for the subsequent analysis. 
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The scalar dominance requires a check of whether the matter contribution term can be ignored at late times. Our 
solution corresponds to a linearly increasing scalar field with positive £ such that the kinetic energy terms and the 
second terms in both (fl~2|) and (fT3| are constant. Therefore, the contribution of the matter density which decreases 
as a -3 ( 1+Wm ) becomes negligible at late times and we ignore the matter part here after. Now, substitution of <j> = t 
leads to 

«W=e* t=f§^\ (A>|). (14) 

In the above equation, we fixed the initial values by a(0) = 1 and 0(0) = 0. Later, we will relax these initial conditions 
and accommodate more general conditions. This describes a de Sitter expansion of the universe. From here on, we 
will always assume a 1 = x 1 and study the time-dependent behavior of the Eqs. (|lip. (|12p. (|13p. Note that for the 
e = +1 case, the cosmological constant term has to be bigger than some positive value. In contrast, for the e = — 1 
case, it could be any value greater than a fixed negative value and we will see that stability adds further restrictions. 
It must especially be negative for stability; still the universe can accelerate and this is driven by the scalar fields. The 
cosmological constant term combines with the contribution from the scalar fields to form an effective cosmological 
constant A eff = 3£ 2 in Eq. (|T4"|) . 



IV. STABILITY 



To check the stability of the above solution, we first consider the following quantities, 

2£<f>-2N = X, N = \n(a) (15) 

Plugging (TTSD into {TTJ) ~ (JT3D , we obtain 

3H 2 + H = -^e x +A (16) 

X + 3HX - (6£ 2 + ^r)e x + 2A = (17) 

The solution (Ti"4"|) corresponds to X — with H = £, 0=1. In order to accommodate more initial conditions, we 
consider the solution X = X(Q) = In/. Then, the effective cosmological constant becomes A e //(/) = — 2e//A 2 + A, 
and the solution (fT4")) is replaced by 



= y/ft + 0(0), a(t) = a(0) e V7« 4 , £ = J-lL + A, (18) 

with 2^0(0) — 21na(0) = In/, A > ^M- . Note that the exponent £ behaves under the change of the initial conditions 

when A = as follows; f —> g, eV^* — * e^*. 
The linear perturbation of Eq. (|T7|) leads to 

SX + 3y/f{6X - (6£ 2 + ^)fSX = (19) 

Introducing SX ~ e' (t , Eq. ([T9j) yields 

7 2 + 3V7^7-(6£ 2 + ^)/ = (20) 

The solutions for the Eq. ([20|) are 



7+/V7 = ^ (21) 



7-/V7 = h (22) 
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Name 


e 




Stability 


A 


(A) 


+1 




7+ > 0, 7- < 0, Unstable, saddle point 


A > 2// A 2 


(B) 


-1 


16/33A 2 < £ 2 < 2/3A 2 


7± < , Stable, attractor 


-6//11A 2 < A < 


(C) 


-1 


£ 2 > 2/3A 2 


7+ > 0, 7_ < 0, Unstable, saddle point 


A > 


(D) 


-1 


£ 2 < 16/33A 2 


7±, Imaginary, stable, attractor 


-2//A 2 < A < -6//11A 2 


(E) 


-1 


£ 2 = 2/3A 2 


7+ = 0, 7- < 


A = 


(F) 


+ 1 


i 2 = 1/3A 2 


Nonperturbative 
4> = V7* + 0o+Aln(l + <7e-V?T' t ) 


A = 3//A 2 



TABLE I: Various accelerating solutions and their stability 



From these equations, we have the following cases: 

(A) . For the e = +1 case, 7 + > and 7_ < for any £. Hence, the solution is unstable. It corresponds to a saddle 
point. The cosmological constant has to be positive. 

For the e = — 1 case, there are three cases. In all three cases, the cosmological constant has to be negative for stability, 
but the late acceleration is achieved by the scalar fields. 

(B) . In the case < £ 2 < j± are both negative and the solution is stable and an attractor. For £ 2 = g^p-, 
the root is degenerate with 7+/V7 = 7-/V7 = — 3£/2. 

(C) . For £ 2 > gjj-, 7 + > and 7_ < 0, so the solution is unstable and corresponds to a saddle point. 

(D) . In the case £ 2 < ggfj, 7± becomes imaginary and the perturbation is oscillatory and it is an attractor. 
The linear perturbation can be integrated explicitly. For (A), (B) and (C), Eq. (fTi?]) yields 

SX = Ae 7 +* + Be 7 "* (23) 

SH = Ael+ * + B "~' - Coe-eV-m), (24) 

A 2 S++6V7£ 7- + 6v7£ ^ K ' 

where Cq = A/(j + + 6\ff£) + -8/(7- + 6\/JQ- Note that when e = +1, the second term in Eq. ([24]) diverges for the 
value f 2 = 1/3A 2 . This might imply that linear perturbation fails in this case and in fact, there exist a nonperturbative 
solution as will be discussed in the next section. For case (D), we have 

5X = L>cos(cj< + o )e"^^ t (25) 



where lu 2 = - 33£ 2 ) and 



where 



6H = -Ee- 6 ^ 1 + F(t)e'iv / /« t (26) 



E = 36L>/V7^cos6l /A 2 (81/^ 2 + 4w 2 ) + 8L>/wsin6» /A 2 (81/^ 2 + 4w 2 ), 
F(t) = 36L>/v / 7?cos(wt + 6' )/A 2 (81/^ 2 + 4w 2 ) + 8Z)/wsmM + 6io)/A 2 (81/e 2 + 4w 2 ). (27) 

We comment on the case £ 2 = 2/3A 2 separately. It corresponds to when the cosmological constant A is zero. In 
this case, one of the roots 7+ of Eq. (f2"Tj) becomes zero, and the other root 7_ is negative. Its stability is indecisive 
at this level. These solutions, their stability and contents of the cosmological constant are summarized in Table I. It 
is interesting to note that in the e = — 1 case, any value greater than a fixed negative of the cosmological constant 
was allowed in (|14l). but stability requires its upper bound must be zero. 



V. NONPERTURBATIVE SOLUTION 



It turns out that in the e = +1 case, an explicit nonperturbative solution can be found. To see that, let us first 
assume = y/fa. Then, Eq. (TP2"]) suggests that for nontrivial solution, we must have 

H 2 = ^\ A 2 A = 3/ (28) 
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The first of the above equation yields £ 2 = 1/3A 2 . This value of £ 2 was the one where linear perturbation failed in 

the previous section. Substituting the ansatz a = -j=e^ into we obtain 

V f 

+ 3£0 2 - 3/^ = (29) 

Note that this equation describes particle motion where constant external force and velocity square dependent frictional 
force are acting. When 0(O) 2 < /, the constant force term dictates the particle motion at early times and it accelerates 
until the velocity reaches the terminal velocity 0(oo) = y/f. When 0(O) 2 > /, the friction term dominates at early 
times and it decelerates until the velocity reaches the terminal velocity 0(oo) = \fj . 
We can find the solution for the above Eq. |29|) as follows 

0(*) = ^t + m-^Hl + C) + ^\n(l + Ce- 6 Vfit), (30) 

with £ = J -Jr? and 



o(t) = a (0)eV«« ( l±g^ffi ». (31) 

The constant C remains arbitrary as long as the validity of the solution is confined within the region 1 + CeV/T* > 0. 
Note that the solution (|30[) has two arbitrary integration constants and it is a complete solution for the ansatz 
= \ffa with 0(0) and 0(0) — v7(l ~ 2C)/(1 + C). It indicates that starting with an arbitrary value of C, the 

solution converges rapidly to cf>(t) = \fjt and a(t) — eV 7 ^*. There is a wide range of initial conditions in which the 
solution rapidly converges to de Sitter acceleration. 

To see the solution more closely, let us divide the case with C > and C < 0. First, note that the condition 

> puts a restriction on the range of |C|, i.e., \C\ < 1 when C < 0. When C > 0, it could be arbitrary. 

We introduce a time scale defined by \C\ — e~ 6 \//£**. Then, Eq. (|3"0|) can be written as 

6(t) = { £ ln (sinh(3V/£(t + <»))) + 0(0), (C < 0, < U < oo) 

| ^ln(cosh(3V7e(i + t»))) +0(0), (C > 0, - oo < i* < oo) 

where 0(0) = 0(0) - ln(sinh(3 % //£t*))/3£ for C < and 0(0) = 0(0) - ln(cosh(3v7£**))/3£ for C > in d32J). And 
we also have 

a(t) = { a(0)(sinh(3V7^ + **))) (C < 0, < t. < oo) (33) 
\ a(0)(cosh(3V7C(* + **))) (C>0, -oo<t*<oo) 

where 5(0) = a(0)/(sinh(3y7^*))^ for C < and 5(0) = a(0)/(cosh(3\/7Ct*))3 for C > in (33]). Note that 0| t= ^ t , 
becomes singular in (|3"2"|) , but our initial time is chosen to be 0, and it is outside the range of dynamics. Had we 
chosen our initial time to be ti, the singularity would occur at t = ti — t*. We will assume that \C\ ~ 1 so that is 
nearly the initial time. 

Let us assume that the initial time t = is chosen when the universe is still at the matter-dominated epoch and 
examine the early-time behavior of the solutions ((32)) and p3|) . Then, for C < with Eqs (|32|) and (|33[) , we have 



<j>{t) - i ln(t + i,) + O , o(t) ~ (t + *.)*. (34) 

The logarithmic time-dependent field [ll|, [l^, EB, also appears in the quintessence with exponential potential. 
The energy density = j?0 2 ~ (t + t*) -2 ~ a -6 scales the same as the stiff matter density and is known as a scaling 
solution. In our case, the scaling behavior holds only at early times, and as time goes by, the full solution (|32p will 
take over. For C > with Eqs. (f32|) and ([33]) . we see that the field and the scale factor a(t) remains constant up 
to first order in time, and as time goes by, both quantities begins to grow. There is no scaling behavior in this case. 
Note that for both cases, Vi^ = — = 1 due to the first condition of Eq. (|2"5)) . This is because in the above solution, 
matter contribution was neglected. It would be interesting to check whether the dominance of the energy density of 
scalar fields emerges from the matter-dominated epoch when the matter contributions are included. 
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Let us discuss some issues of dark energy with this solution. First, we have the acceleration given by 

a ^ „l±10 e - 6 V/«( t+t -) +e- 12 V7« t +*.) 

- = /£ 7= , (35) 

a (i ± e -6v7S(t+*,))2 

where + in the numerator is for C > and — is for C < 0. We see that for C > 0, it is always accelerating. For 
C < 0, there is a transition time where the scale factor changes from deceleration to acceleration. It is given with 
i* ~ by 

^ = -±=- ln(5 + 2>/6) ~ 0.48i . (36) 

where to is the current age of the universe [32]. This fix v7C ~ 0.79/io- Next, equation of state with this value of 
V/£ is given by 

e -6V7«(*+*«) 

cj = -1=f8 = It 0.07. (37) 

(1 zp e - 6 V/«(*+**))2 

For C > with an upper sign in (|37[). cj approaches —1 from below, w is singular at i = — and again it is outside 
the range of dynamics. For C < 0, it approaches —1 from above. And there is no singularity. In both cases, there is 
no crossing the w = — 1, but asymptotically approaches the uj = — 1 line. We find that this number for the equation 
of state in (|37|) is in the region allowed by the observational data Q , but inclusion of matter might change the result 
somewhat. 



VI. CONCLUSION AND DISCUSSION 



We presented a de Sitter nonlinear sigma model coupled to Einstein gravity in order to describe the current 
acceleration of the Universe. It has some characteristic features as follows. Out of the four scalar fields, only one of 
them is time-dependent and the remaining three fields have only spatial dependence. If the time dependent scalar field 
is phantom, then the remaining fields are ordinary scalar fields, or vice versa. The formal case could also be thought 
of as the dilatonic phantom coupled with triplet of scalar fields, whereas the latter case as the dilaton coupled with 
triplet of phantom scalar fields. Since kinetic energy of both positive and negative sign exists, it could be thought 
of as a quintom model [268 ] with dilaton interaction between the two sectors. But the quintom model only considers 
time-dependent fields. A specific form of the potential is not needed to achieve the late-time exponential acceleration, 
but introduction of a potential could produce subdominant behavior of the scalar fields as in the quintessence. It is 
suspected that the potential does not modify the late-time exponential behavior, because it does not change the Eq. 

We find that a simple ansatz provides the constant energy density and results in an accelerating universe with 
an exponential expansion. The balance between the pressures coming from the time-dependent field and spatial- 
dependent fields makes it possible to achieve the exponential acceleration. It is pointed out that the target space of 
Euclidean de Sitter space with signature (+, +, +, +) cannot produce such balance and exponential acceleration of the 
universe. The model has essentially two parameters, £ and the cosmological constant term A. y/j£M p plays the role of 
the Hubble constant and is a function of the strength of the self-coupling constant and the cosmological constant term 
A. Consider, for example, the nonperturbative case with £ 2 = 1/3A 2 and A = 3//A 2 . Recall / = e 2 €0(o)/M p / fl 2(o) 
and let us assume a(0) ~ 1. We mention a couple of cases where y/j£M p ~ 10 _61 M p and A ~ 10~ 122 M 2 can be 
realized. In the first case with (f)(0) ~ — M p and A ~ 1/250, we have £4>(0) ~ — H5M p . If the scale when the nonlinear 
sigma model sets in is of the order of Gev with cj>(0) ~ —lGeV, this requires extremely weak coupling constant with 
A ~ 10~ . In this case, we have £</>(0) ~ — 190M P . These are fine tunings which can yield the small Hubble constant 
and the cosmological constant. 

Stability analysis shows some of the solutions, depending on the values of the parameter £, are classically stable and 
attractor solutions. They require that the original cosmological constant term must be negative, still the acceleration 
is possible led by the scalar fields. In one case, where the cosmological constant term is uniquely fixed, there is 
a nonperturbative solution which asymptotically approaches the de Sitter phase of acceleration. This solution also 
exhibits a power law expansion at early times, and the energy density of the scalar fields mimics the matter energy 
density. It remains to be seen whether the stability survives when the analysis is extended to spatial variations. 

The present analysis indicates that the acceleration phase can be dominated by the nonlinear sigma model. We 
only focused on the late time behavior except the nonperturbative case. To show whether this behavior of scalar 
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dominance can emerge from matter-dominated epoch, the analysis has to be extended including the contribution 
of matter density at early times which was neglected. Finally, whether the de Sitter nonlinear sigma model could 
come from particle physics as an effective low energy field theory remains to be seen. These aspects needs further 
investigation. 

Note added: After the completion of this work, we became aware of Ref. [33| where the ansatz ([5]) each multiplied 
by some constant factors to have Minkowski background also appeared in the cosmological context of the Lorentz 
violating massive graviton models [33] . These models deal with flat background metric. However, in our de Sitter 
background solution, a linear perturbation of the metric in the action ([1]) does not result in any massive graviton 
mode even though the ansatz ([5]) spontaneously break the diffcomorphism invariance. This can be readily seen by 
checking that the mass term which is of the second order in the perturbations cancels out in the second term of the 
action in ([T]). Moreover, the de Sitter background solution allows the modification where each of the ansatz in ([5]) can 
be multiplied by a same constant factor only which still does not yield Lorentz violating mass term. Perhaps, it could 
be possible to generate a mass term by a suitable deformation of the target space metric in the action {!]), which is 
beyond the scope of this work, but nevertheless whose implications would be worthwhile to be explored in detail. 
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